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We analyze the stability of excitonic ground states in the two-band Hubbard model with addi¬ 
tional electron-phonon and Hund’s rule couplings using a combination of mean-field and variational 
cluster approaches. We show that both the interband Coulomb interaction and the electron-phonon 
interaction will cooperatively stabilize a charge density wave (CDW) state which typifies an “ex¬ 
citonic” CDW if predominantly triggered by the effective interorbital electron-hole attraction or a 
“phononic” CDW if mostly caused by the coupling to the lattice degrees of freedom. By contrast, the 
Hund’s rule coupling promotes an excitonic spin density wave. We determine the transition between 
excitonic charge and spin density waves and comment on a fixation of the phase of the excitonic 
order parameter that would prevent the formation of a superfluid condensate of excitons. The im¬ 
plications for exciton condensation in several material classes with strongly correlated electrons are 
discussed. 

PACS numbers: 71.10.Fd, 71.35.-y, 71.45.Lr, 71.30.-l-h 


I. INTRODUCTION 

Excitonic insulator (El) phases show a spontaneous co¬ 
herence between conduction-band electrons and valence- 
band holes, where the prior formation of bound electron- 
hole pairs (excitons) is typically triggered by the inter¬ 
band Coulomb interaction. The condensation of excitons 
in an El state was theoretically predicted half a century 
ago to occur in semiconductors (semimetals) with small 
band gap (band overlapThe condensed excitonic 
phase can be characterized by a nonvanishing order pa¬ 
rameter (cj[,_|_Q/fc), where cj[ and fl are the creation oper¬ 
ators of an electron in the conduction and valence bands, 
respectively. If the valence-band top and conduction- 
band bottom are separated by the wave vector Q, the 
system forms a density wave with modulation Q. It 
is important to note that the El—despite representing 
a macroscopic, phase-coherent quantum state—does not 
necessarily feature supertransport properties. 

The experimental efforts to establish the El in weakly 
correlated bulk materials largely failed. It is only re¬ 
cently that exciton condensation has been addressed in 
systems with rather strong electronic correlational. In 
this regard, Tm(Se,Te) was argued to exhibit a pressure- 
induced excitonic instability, related to an anomalous in- 
cr ease in the electrical resistivity and thermal diffusiv- 
itjian. The charge-density-wave (CDW) state observed 
in lT-TiSe2 was claimed to be of excitonic origirPISl 
Cai_a;Laa;B 6 , the weak ferromagnetism was interpreted 
in terms of doped spin-triplet excitonJi^Hi^. The con¬ 
densation of spin-triplet excitons was also predicted to 
occur in the proximity of the spin-st ate tr ansition,!^ of 
which Pro.sCao.sCoOa is an example.l^^Iini Likewise, the 
structural phase transition of the layered chalcogenide 
Ta2NiSe5 has been attributed to a spin-singlet El 21 Il 2 ^. 
The spin-density-wave (SDW) state of iron-pnictide su¬ 
perconductors has sometimes been argued to be of the 


excitonic origin as welP^^^. Einally, an El state was 
suggested in a t2c)-orbital system with strong spin-orbit 
couplinj^. 

From the theoretical side, the extended Falicov- 
Kimball model was considered as a paradigmatic model 
to describe the El formation and the closely related phe¬ 
nomenon of electronic ferroelectricitjl 22 MlJ Here the spin 
degrees of freedom were not taken into account, how¬ 
ever. Excitonic phases in strongly correlated spinful sys¬ 
tems can be discussed in the framework of two-band 
Hubbard-type models. Including thereby the Hund’s rule 
coupling is known to stabilize the spin-triplet excitonic 
phase in the otherwise de generate s pin-singlet and spin- 
triplet excitonic phaseji^l^lllllllll! Qn the other hand, 
we have shown in our previous worlP^ that taking into 
account electronic interactions only, a spin-singlet ex¬ 
citonic phase cannot be stabilized, which may however 
be realized in lT-TiSe2 and Ta2NiSe5, where the impor¬ 
tance ofdectron-phonon coupling was recently pointed 
oulfi^HIillSl. Although the spin-singlet excitonic state 
has been investigated in the spinles s multiband model 
with electron-phonon couplin^i^llllll]^ not much is known 
about the role of the electron-phonon coupling played in 
the excitonic density wave states in the spinful multiband 
Hubbard model. 

Motivated by the recent developments in the held, in 
this paper, we will thoroughly investigate the stability of 
the excitonic density wave states in the two-band Hub¬ 
bard model with additional electron-phonon coupling and 
Hund’s rule exchange. The model is analyzed employing 
static mean-held theory for the electron-phonon coupling 
and the variational cluster approximation (VGA) for the 
electronic correlations. In doing so, we will hrst show 
that the interband Coulomb interaction U' and electron- 
phonon interaction A cooperatively stabilize the CDW 
and that a smooth crossover occurs between “excitonic” 
CDW and “phononic” CDW states, just by increasing 
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the ratio XjU'. Incorporating the Hund’s rule coupling 
J, an excitonic SDW state competes with the excitonic 
CDW. The ground-state phase diagram of such extended 
two-band Hubbard model is determined in the J-A plane. 
We will, moreover, pay particular attention to the phase 
of the order parameter in the presence of the electron- 
phonon and Hund’s rule couplings and show that both 
electron-phonon coupling and pair-hopping terms fix the 
phase of the excitonic order parameters, thereby prevent¬ 
ing the system from realizing a superfluid. Finally, the 
implications for exciton condensation in real materials 
will be discussed. 

The paper is organized as follows: In Sec. H, we present 
our model and briefly outline the methods of calculations. 
The numerical results will be presented in Sec. HI, where 
electron-phonon interaction, Hund’s rule coupling, and 
pair-hopping effects will be discussed and the ground- 
state phase diagram is derived. Section IV relates our 
results to recent experiments and draws conclusions. 


II. MODEL AND METHODS 

A. Model Hamiltonian 

We consider the two-band Hubbard model, supple¬ 
mented by electron-phonon and Hund’s rule coupling 
terms, 

n = ne + nl, + ni, + nph + Ue.ph , (i) 

defined on a two-dimensional square lattice. The nonin¬ 
teracting //c-band electrons are described by 

^fc/fecr/fco- + (2) 

h.<7 k,(7 

where {^ka) denotes the creation operator of an 
electron with momentum k and spin a (=t,4,) in the 
/ (c) band. Within the tight-binding approximation, 
the dispersion of band a (= /, c) is given by = 
—2ta{coskx + cosky) + Ea — fJ-, where ta is the electron 
hopping integral between the neighboring sites and Ea is 
the on-site energy of the a orbital. We assume £/ < 0 
and £c > Oj so that the / and c bands correspond to the 
valence and conduction bands, respectively. The chem¬ 
ical potential ^ is fixed to ensure a filling of two elec¬ 
trons per site (half filling), i.e., (n{) + {nf) = 2 with 

< = 

The repulsive Coulomb interaction takes the form 

i i i 

where Ua is its intraorbital part and U' gives the in¬ 
terorbital contribution that is responsible for an effective 
electron-hole attraction and eventually for an excitonic 


instability in the system. The Hund’s exchange interac¬ 
tion is defined by 

nU=-2jY^{Sl-S<l + \n{n!l) 

i 

-J'Y. + 4t4ifii-fii) ( 4 ) 

i 

with S°‘ = a' ot'i^cTaaiOLia' /2, where <J is the vector of 
Pauli matrices. J and J' are the strengths of the Hund’s 
rule coupling and pair-hopping term, respectively. J and 
J' will stabilize a spin-triplet excitonic stated!] 

In Eq. 0 , we also included the phonon degrees of free¬ 
dom because the lattice displacements play an important 
role in the materials under consideration. The electron- 
phonon coupling becomes particularly important when 
we address spin-singlet electron-hole excitations. In the 
harmonic approximation, the phonon part of the Hamil¬ 
tonian is given by 

— 'y ) ( 5 ) 

q 

where the bosonic operator fol creates a phonon with mo¬ 
mentum q and frequency coq (we have set H = 1). The 
dominant electron-phonon coupling term between a c-f 
(electron-hole) excitation and lattice displacement is as¬ 
sumed to be 

He-ph = —^ + H-C. (6) 

^ cr 

with coupling constant 

Throughout the paper, we fix the hopping parameters 
tf = tc = t and use t as the unit of energy. Furthermore, 
we set Ec/t = —Efjt = 3.2, so that the noninteracting 
band structure represents a semimetal with a small band 
overlap. The conduction-band bottom ai k — (0,0) gives 
rise to an electron pocket, while the valence-band top 
produces a hole pocket at A; = (tTjTt), resulting in the 
modulation vector of the density wave Q = (tt, tt); see 
Ref. 02] for the band dispersion and Fermi surface in the 
Brillouin zone of the square lattice. For simplicity, we 
assume Uf = Uc = U and employ U = 2U' — J to sup¬ 
press the Hartree shiff!^. In this choice, the FI state is 
stabil ized b etween the band-insulator and Mott-insulator 
state^42143l_ Moreover, we consider a dispersionless Ein¬ 
stein phonon Wg = w and a momentum-independent 
electron-phonon coupling constant gq = g. Since the 
strength of the electron-phonon coupling appears in the 
form A = 5 ^/a; in the mean-field approximation used be¬ 
low, we take A as the electron-phonon coupling parameter 
in what follows. 


B. Mean-fleld approximation for the phonons 

We treat the electron-phonon interaction term He-ph 
in the mean-field (frozen-phonon) approximation. In- 
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troducing the expectation values of the c-f hybridiza¬ 
tion (cV) and lattice displacement (6), the opera¬ 
tors in Eq. ^ are approximated as ~ 

[ (^q)cI+q(j/fcCT+^q(Cfc^qo'/fco') ]^q,Q~{^q) {^k+qa-fko-)^q,Q- 

Since in our model the nesting vector Q = (tt, tt) is com¬ 
mensurate with the lattice periodicity, = 1 for lat¬ 

tice vectors r^. Hence, we have bq = b_Q {bq = 
where bq and b-q (bq and b^q) annihilate (create) 
the same phonon. This implies (bq) = (b-q) = (bq), 
and therefore (bq) becomes a real number. In view of 

(^i+Qafk<y) = ifla^^k+Qa)* 7 ^ 0, we define the complex 
order parameter of the excitonic CDW as 

= = (7) 

k.CT 

where |<i>c| and 9c are the amplitude and phase of the 
order parameter, respectively. Then the electron-phonon 
part in the mean-field approximation is 


expectation values have opposite signs on the same 
site. In momentum space, this reads X)fe('^I+Q'f/fcT) = 
“ I]fc(cI+Q;/fcf)- We then find, from Eqs. 0 and 

that (bq) = (6^) oc Efc,^(4+Qa/fe^) = which means 
that the spin-triplet condensate will not couple to the 
phonons. 

C. Variational cluster approximation 

In order to take electron correlation effects into ac¬ 
count, we treat the electronic interactions in Q within 
the VCA,E3311 which is a quantum cluster method based 
on the self-energy functional theorjB^. The VGA first 
introduces disconnected clusters of finite size, for which 
the cluster self-energy E' can be computed exactly. In a 
next step, out of this, a superlattice is formed as a refer¬ 
ence system. By restricting the trial self-energy to E', we 
obtain a certain approximation to the grand potential of 
the original system. 


'H^ph = X! i/k+Qafka + fLck+qa^ 

k.a 

+ -f bq)\^c\cos 9c - 8gVN{bq)\^c\ cos 0c. 

( 8 ) 


Introducing Bq = bq + Sq^q(4:gVN/u})\^c\cos9c, the 
phonon Hamiltonian Bph together with the second term 
of the right-hand side of Eq. ® can be diagonalized, 
yielding ojY.qB\Bq - 16AV|$^ cos^ 9c- Hence, from 
{Bq) = (H^) = 0, we find 


{bq) = {h^q)=-^^mcos9c. (9) 

Substituting this expression into Eq. (§, we finally ob¬ 
tain the mean-field electron-phonon Hamiltonian, 


cos 9c ^ fLck+qa + H.c. -f 

k,(T 


NAl 

4A 


cos^ 9c 


( 10 ) 


with Ap = —8A|<I>c|. Using Eq. (10), below we will min¬ 
imize the grand potential of the system with respect to 
Ap and 9c- 


We define the complex order parameter of the excitonic 
SDW as 




1$. 


Ue. = 


1 


E ^(Cfc+Qo-ZfcCT), 

k,(7 


( 11 ) 


where |4)s| and 9s are the amplitude and phase of the 
order parameter, respectively. Because we assume an 
SDW state with modulation vector Q — (tt, tt), where 
the expectation value (c^^/i^) is in antiphase compared 
to regarding the spatial variation, these two 


+ Ti ln(Go 1 - E')”^ - Tr In(G'), (12) 

where and G' are the grand potential and Green’s 
function of the reference system, respectively, and Go is 
the noninteracting Green’s function; for further details, 
see Refs. 1501511 In doing so, the short-range electron 
correlations within the cluster of the reference system 
are treated exactly. In our VGA calculation, we take an 
Lc = 2x2 = A site (eight-orbital) cluster as the reference 
system and we use exact diagonalization to solve the cor¬ 
responding quantum many-body problem in the cluster. 
Within VGA, we can take into account spontaneous sym¬ 
metry breakings just by adding appropriate Weiss fields 
to the reference systenPH^ and take these fields as vari¬ 
ational parameters. The Weiss fields for excitonic CDW 
and SDW states, which are defined by the order parame¬ 
ter [in Eq. 0 ] and 4>s [in Eq. respectively, may 

be written as 

Hr = A'e*®=^/Lcfc+Q. + H.c. (13) 

k,a 

= Ay®=^a/Lcfc+Q.+H.c.. (14) 

k,(7 


Here, Ag and A), are the strengths of the Weiss fields for 
the excitonic GDW and SDW states generated by 
and 


According to Eq. (10) , we take into account the con¬ 


tribution of the phonons in the mean-field approximation 
as a one-particle term in the original system. Then, the 
Hamiltonian describing an excitonic GDW state in the 
reference system is given by 


n) = He+ n 


u 

e-e 


. I '^MF 

^e-e ' ^e-ph 


K 


WF 


(15) 


where we note that 'He+'Hc_c+'H'c_c + 'H^pi^ is the Hamil¬ 
tonian of the original system and the Weiss field 
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is added in the reference system. Using we calcu¬ 
late the grand potential U and optimize the variational 
parameters Ag, Ap, and 0c- The most stable solution 
with (Ag, Ap) ^ (0,0) corresponds to the excitonic CDW 
state. Note that we determine the parameters Ap and 9c 
via the minimization of the grand potential rather than 
solving the self-consistent equation. Both procedures are 
equivalent, however, since the order parameter <i>c cal¬ 
culated, using the Green’s function with Ap and 0c opti¬ 
mized via the grand potential calculation in VGA, exactly 
satisfies the self-consistent condition Ap = 8A$c. 

Since the spin-triplet term does not couple to the lat¬ 
tice degrees of freedom within our mean-field approach, 
the phonons will not affect the excitonic SDW state. 
Then the Hamiltonian of the reference system describ¬ 
ing an excitonic SDW is 

K^He + Hlc + nic + (16) 

Again we calculate the grand potential H from the ref¬ 
erence Hamiltonian T-L'^ and optimize variational param¬ 
eters Ag and 0s, where the most stable solution with 
Ag 7 ^ 0 corresponds to the SDW state. 


III. NUMERICAL RESULTS 


A. Phase of the order parameters 


We first discuss the phase of the different order param¬ 
eters entering the grand potential. In the spin-singlet ex¬ 
citonic state, the system forms an excitonic CDW at any 
finite U' and A due to the perfect nesting of the Fermi 
surface. Figure [^a) shows the calculated grand poten¬ 
tial U as a function of the variational parameters Ag and 
Ap. Obviously the grand potential has a stationary point 
at (Aq,Ap) ^ (0,0), signaling a CDW ordering. With¬ 
out electron-phonon coupling, fl is independent of the 
phase 0c, i.e., U(0c) = ^{0c)- Accordingly, the excitonic 
CDW state reveals a^apless acoustic phase mode in its 
excitation spectruni^. If, however, the electron-phonon 
coupling comes into play, the grand potential manifests 
a dependence on the phase of the (complex) order pa¬ 
rameter. In Fig. [^b), we display the 0c dependence of 
U; the grand potential takes its minimum at 9c = 0,7r. 
This phase fixation may be expected looking at Eq. (101. 


In our mean-field approximation, the single-particle gap 
caused by A is given as Ap cos 0c and is maximized at 
0c = 0. When 0c is fixed by the electron-phonon cou¬ 
pling, the collective phase mode in the spin-singlet ex¬ 
citonic state becomes massive (see the discussion of the 
spinless model in Ref. l45ll . 

In the case of the spin-triplet excitonic state, the 
excitonic SDW and CDW states are degenerate if the 
electron-phonon and Hund’s couplings are neglected. 
The Hund’s exchange terms oc J and oc J' lift this de¬ 
generacy and stabilize the SDW stat^^. Note that the 
9s dependence of the grand potential behaves differently 




FIG. 1: (Color online) (a) Grand potential H as a function 
of the variational parameters Ag and Ap for U'/t = 4 and 
\/t = 0.15. (b) 9c dependence of H taken the values of Ag and 
Ap optimized at 9c = 0. Green dots indicate the stationary 
points. 



FIG. 2: (Color online) (a) Grand potential H as a function of 
the variational parameter A( at U'/t = 4 and Jjt = J'/t = 
1. (b) 9s dependence of Q. obtained using the value of A), 
optimized at 9s = 0. Dots mark stable stationary points. 


in the presence or absence of the pair-hopping term J': 
For J' = 0, the grand potential of the SDW state does 
not depend on 0s, i.e., H(0s) = 11(0'), whereas fl depends 
on 0s at any finite J'. Again the independence of 4>s on 
the phase value 9s accounts for a gapless excitation spec¬ 
trum, i.e., an acoustic phase mode. Figure gives the 
calculated grand potential fl as a function of the phase 
0s in the presence of the pair-hopping term J'. Indeed 
we find that H has two minima, at 9s = 0, tt, which fixes 
the phase 0s of $g. It is known that the energy in the 
presence of the pair-hopping-type exchange interaction 
shows a phase dependence cos20j^. This is why the 
pair-hopping term J' fixes Og and in that way destroys 
the gapless acoustic phase mode in the spin-triplet exci¬ 
tonic state. 


B. Excitonic CDW state 

Now let us analyze the stability of the CDW state in 
the presence of the electron-phonon coupling in more de¬ 
tail. In Fig. 1^ we present the results for both the op¬ 
timized grand potential Hopt and the order parameter 
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FIG. 3: (Color online) Optimized values of the grand potential 
flopt in dependence on (a) U' jt and (b) A/t. Here, Hq is the 
grand potential in the normal (semimetallic) state. Order 
parameter for the excitonic CDW state as a function of 
(c) \J'jt and (d) Xjt. 


when the interband Coulomb interaction U' and the 
electron-phonon coupling A are varied, flopt indicates 
that (i) the symmetry-broken CDW state is lower in en¬ 
ergy than the normal state and (ii) the stability of the 
CDW state is enhanced if U' and A are increased; see 
Figs. ia) and[^b). This is corroborated by the behav¬ 
ior of the order parameter displayed in Figs. [^c) and 
id). We see that the interband Coulomb interaction U' 
induces and boosts the excitonic CDW state while the 
electron-phonon coupling A rather promotes a phononic 
CDW state (see below). Both, however, cooperatively 
stabilize a charge-ordered state. In this connection, the 
electron-phonon coupling lifts the degeneracy of CDW 
and SDW that exists for A = 0. 

In the mean-field approximation, the gap parameter of 
the CDW state, Ac = (17'-|-8A)$c, can be separated into 
two contributions: the excitonic (or interband Coulomb 
driven) part Aq = U'^c and the phononic (or electron- 
phonon driven) part Ap = 8A$c. Figure ®a) illustrates 
the relative magnitude of Aq and Ap, in dependence on 
the ratio 8X/U'. At 8X/U' ^ 1, Ac ~ Aq ^ Ap and the 
CDW state, stabilized by the interband Coulomb inter¬ 
action U', is excitonic by its nature. Increasing 8X/U', 
Aq decreases while Ap increases, indicating a smooth 
crossover to a phononic CDW, which fully develops at 
8X/U' I, where Ac ~ Ap Aq. In the crossover 
region 8X/U' ~ I, both excitonic and phononic contribu¬ 
tions are equally important. 

In Figs. [^b)-[^e), we show the behavior of the dif¬ 
ferent contributions to the gap parameter Ac when U' 



FIG. 4: (Golor online) (a) Excitonic gap parameter Aq (solid 
line) and phononic gap parameter Ap (dashed line) divided by 
the total gap Ac = Aq -I- Ap. Ac, Aq, and Ap are separately 
plotted as a function of (b),(c) U'/t and (d),(e) X/t. 


and A are varied separately. Data are obtained by VCA. 
Enhancing U'jt {Xjt) at weak Xjt (small Vjt') leads 
to an increase in Ap (Aq) as well, since both interac¬ 
tions couple to the same operator-product expectation 
value (cj[^Qcr/fccr); see Figs, dlb) and 4ld). The crossover 
between excitonic and phononic CETws can be seen in 
Figs.j^c) and|^e), where a crossing between Ap and Aq 
appears whenG' ~ 8A. 


C. Excitonic SDW state 

We now study the influence of the Hund’s rule cou¬ 
pling on the nature of the excitonic phase, and also when 
an additional electron-phonon coupling acts in the sys¬ 
tem. Evidently, excitonic CDW and SDW states are de¬ 
generate at J = J' = 0 and A = (PHHl Any finite J 
and/or A lifts this degeneracy. Figure clearly shows 
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FIG. 5: (Color online) J dependence of the grand potential 
flopt and the order parameter "F in the excitonic CDW (sym¬ 
bols) and SDW (solid line) states (a),(c) with {J' = J) and 
(b),(d) without (J' = 0) the pair-hopping term, where U'/t = 
2.4. flo is the grand potential of the normal semimetallic 
state. 



FIG. 6: (Color online) Ground-state phase diagram of the 
two-band Hubbard model with electron-phonon and Hund’s 
rule couplings showing the stability regions of excitonic CDW 
and SDW phases. Results are obtained, in the presence 
(J' = J) and absence {J' = 0) of the pair-hopping term, 
by combining the mean-held and VGA approaches, for a two- 
dimensional (square) lattice at half Filing, where U'/t = 2.4. 


that by increasing J, the optimized grand potential ilopt 
for the SDW (CDW) state monotonically decreases (in¬ 
creases); accordingly, the order parameter for the SDW 
(CDW) phase is enhanced (suppressed). This holds for 
both J' > 0 and J' = 0. Clearly the SDW state is sta¬ 
ble as soon as becomes less than A finite 


pair-hopping term oc J' amplifies the tendency towards 
SDW formatioiP^. 


D. Ground-state phase diagram 

The competition between electron-phonon and Hund’s 
rule coupling effects leads to the ground-state phase di¬ 
agram of the model 0 presented in Fig. Obviously, 
A and J tend to establish CDW and SDW phases, re¬ 
spectively, on top of an excitonic state enforced by U'. 
A finite J' increases the region in the J-A plane where 
the excitonic SDW is the ground state. We note that the 
SDW-CDW transition is a first-order transition, within 
the limits of our approximations. 


IV. DISCUSSION AND CONCLUSIONS 

First, let us discuss implications of our findings on 
materials aspects. The transition-metal chalcogenides 
ir-TiSe 2 and Ta 2 NiSe 5 have recently been discussed in 
terms of the spin-singlet EL In these systems, the valence 
and conduction bands are formed by orbitals located on 
different atoms. For example, in lT-TiSe 2 , the Ap or¬ 
bitals of Se ions account for the valence bands and the id 
orbitals of Ti ions account for the conduction band^^®^, 
and in Ta 2 NiSe 5 , the id orbitals of Ni ions form the 
valence bands and the id orbitals of Ta ions form the 
conduction band^^Ut^ll Hund’s rule coupling, acting be¬ 
tween electrons on different orbitals of a single ion and 
favoring the spin-triplet excitons, is therefore negligible. 
Rather, in these materials, the electron-phonon coupling 
is at play and will stabilize a spin-singlet El state. The 
interband Coulomb interaction and electron-phonon in¬ 
teraction, which are inherently interrelated in these ma¬ 
terials, will cooperatively stabilize the El CDW, which 
is predominantly phononic or excitonic depending on the 
importance of electron-phonon or Coulomb effects. 

By contrast, in the iron-pnictide superconductor^^SMl 

and Co oxideJi^H^, the valence and conduction bands 
are formed by the d orbitals on the (same) transition- 
metal ions, so that the Hund’s rule coupling is expected 
to be strong. Hence, in these materials, the SDW phase, 
if really excitonic in origin, is rather triggered by the 
Hund’s rule coupling than by electron-phonon coupling. 
Then, as our phase diagram suggests, the condensation 
of spin-triplet excitons will play a major role. 

Second, let us comment on the phase of the exci¬ 
tonic order parameters. On the one hand, as we have 
shown in the preceding section, the electron-phonon in¬ 
teraction stabilizes the spin-singlet excitonic condensate, 
whereas exchange interactions such as the Hund’s rule 
couplings stabilize a spin-triplet excitonic condensate in 
the otherwise degenerate excitonic density-wave states. 
On the other hand, these interactions, in particular the 
electron-phonon and pair-hopping interactions, will fix 
the phase of the order parameter of the excitonic state; 
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see Sec. Ill A. Because the spatial modulations of the 
CDW and SDW are given by cos(Q • + 0), the phase 9 

may lead to a translational motion of the condensate as 
a whol^^. If the energy of the condensate is independent 
of the phase, maintaining the continuous symmetry of the 
system with respect to the phase, a gapless acoustic phase 
mode may appear in the excitation spectrum, allowing 
for a translational motion of the condensate without loss 
of energy (i.e., superfluidity), as predicted by Frohlich in 
his theory of incommensurate density wave^^. In real 
materials, however, excitonic condensation will be influ¬ 
enced by the lattice degrees of freedom or affected by the 
pair-hopping term. Then the phase of the condensate 
is fixed and a gap opens for the collective phase mode. 
This makes realization of excitonic superfluidity in real 
materials unlikely. 

To summarize, we have studied the stability of the ex¬ 
citonic states with charge and spin density modulations 
in terms of the two-band Hubbard model, supplemented 
by electron-phonon and Hund’s rule interactions, where 
the static mean-field theory is employed for coupling to 
the lattice degrees of freedom and the variational cluster 
approximations for the electron correlations. We have 
shown that both the interband Coulomb interaction U' 
and the electron-phonon coupling A tend to stabilize an 
excitonic CDW state. While at A = 0 the excitonic insu¬ 


lator exhibits an acoustic phase mode, any finite A fixes 
the phase of the order parameter and therefore eliminates 
such a gapless excitation related to supertransport prop¬ 
erties. The CDW typifies a predominantly excitonic and 
phononic state for small and large ratios 8X/U', respec¬ 
tively. The Hund’s rule coupling J, on the other hand, 
promotes an excitonic SDW phase, which is further sta¬ 
bilized by pair-hopping processes, which also fixes the 
phase of the order parameter. These results obtained 
for a generic microscopic model Hamiltonian should con¬ 
tribute to a better understanding of exciton condensation 
in several material classes with strong electronic correla¬ 
tions. 
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